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ABSTRACT 


~”= 


A two-dimensional finite eiement model of a low frequency, 
free-flooded piezoelectric thin ring sonar transducer was 
developed for use with the ATILA code. 

Effective material properties for the ring were determined 
from in-air modal analyses using ATILA. These were adjusted to 
obtain the closest agreement between the calculated ring 
resonance frequency and coupling coefficient and measured data 
supplied by the manufacturer. In addition, equations for the 
electrical admittance for this transducer in air were 
developed based on an analytical approach published by Hong- 
Zhang Wang (J. Acoust. Soc. Am. 79, 164-176 (1986)). A Matlab 
program was written to implement these equations and to 
calculate various electromechanical network parameters. The 
calculated values agreed with measured values to within 5 
percent. 

An in-water harmonic analysis, in which the transmitting 
voltage response and the directivity pattern are computed, was 
performed using the ATILA code. In general, results obtained 
using the finite element model and manufacturer’s measured 
values agree within 10 percent. The effects of changes in 
material properties on acoustic performance of the transducer 


were investigated; the results are discussed. 
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I. INTRODUCTION 


rree-flooded piezoelectric ring transducers have found 
increasing use in deep ocean sonar applications. Thougn simpie 
in concept, the ability to reliably predict the performance of 
fem. Gevicses iS Complicated by the nature of their 
@emstruction. In particular, the acoustic performance of such 
Beemcducers depends upon the physical properties of the 
encapsulant materials used, making their performance difficult 
ce model. 

Standard multiport network analysis techniques, long used 
to model the electromechanical behavior of piezoelectric 
transducers, are usually restricted to one-dimensional motion. 
Newer types of transducer structures, such as the free-flooded 
ring, that have more than one important degree of freedom, do 
not lend themselves to this technique. Consequently, finite 
element analysis (FEA) techniques, which provide’ the 
capability to analyze engineering problems of many degrees of 
Freedom, are frequently employed. 

The objectives of this research were (1) to numerically 
model a free-flooded piezoelectric ring sonar transducer and 


(2) using the developed model, predict the changes in acoustic 


perrormance due to changes in the elastic properties Giga. 
encapsulants used. 

The electroacoustic properties of this free-rileeeea 
piezoelectric ring sonar transducer were numericaily modeled 
uSing the finite-element code "ATILA", developed by the 
Institut Supérieur d’Electronique du Nord (Lille, France) 
[Ref. 1]. Since the ring is constructed of ceramic staves 
cemented together, the effective material properties for the 
ring were determined from in-air modal analyses using ATILA. 
The material property values were adjusted to obtain the 
closest agreement between the calculated ring assembly 
resonance frequency and coupling coefficient and measiredsaaws 
Supplied by the manufacturer. 

In addition to the finite element analysis, equations for 
the electrical admittance of the thin ring transducer when 
unloaded were developed based on an analytical approach 
published by Hong-zhang Wang [Ref. 2]. A MATLAB program was 
written to implement and compute various electromechanical 
network parameters. The results obtained from the analytical 
method were compared with the measured values and those 
obtained by the finite element code. 

The finite element model built was further used to predict 
the changes in acoustical properties due to a change in the 
type of encapsulant materials used in manufacturing. 

The remainder of this thesis is divided into seven 


chapters. Chapter II describes the transducer under study. 


Chapter III gives a generai description of the finite element 
analysis and the standard network analysis of piezoelectric 
ring transducers. Chapter IV describes the development of tne 
2-dimensional finite element model of the transducer tor use 
with the ATILA code. Chapter V presents and discusses the 
results of the in-air modal analyses and in-water harmonic 
analyses using ATILA, as well as the results of the analytical 
model. Chapter VI presents the conclusions. Chapter VII 
provides suggestions for future work. Appendix A contains a 
copy of the input data file for the finite element mesh used 
in the in-water harmonic analyses. The MATLAB program written 
to implement the analytical model for the transducer in air is 
attached as Appendix B. Appendix C provides a list of the 
typical values of material constants used in the modeling. 
Appendix D is a sample calculation of the effective elastic 
constants of a homogeneous ceramic ring incorporating the 


fiberglass wrap. 


Il. TRANSDUCER DESCRIPTION 


Figure 1 shows a schematic diagram of a free-flooded 
plezoelectric ring sonar transducer. When immersed in water 
and vibrating radially, such a transducer exhibits two 
fundamental resonant modes that can radiate acoustically into 
the water. They are the radial resonance of the ring and the 
Cavity resonance of the enclosed water column. Radial motion 
of the cylinder walls will excite the symmetrical cavity modes 
in the enclosed water column. The radial resonance is the more 


important resonance. 





Figure 1: The Free Flooded Ring 
Transducer 


The moror element (piezoelectric ring) of the transducer 
under study consists of a set cf 72 tangentially volarized 
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1@aq zirconate titanate, PZT, .Navy Type I} ceramic staves as 
shown in Figure 2. 

PZT 1S an isotropic poliycrystaiine ceramic which exhibits 
piezoelectric properties after poling [Ref. 4]. It exhibits a 
nign degree of symmetry in the plane transverse to the poling. 


These materials are usuaily treated as member of the 6 mm 


hexagonal class [Ref. 4]. 
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Figure 2: The Ceramic Staves 


Each stave has dimensions (wxl1xh) 0.25 x 0.6519 x 6.8 
in. (0.635 x 1.656 x 17.272 cm). They are cemented together 
using epoxy resin to form a ring of 14.44 inches (0.3668 m) 


outside diameter, as shown in Figure 3. 
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Figure 3: The Ring Transducer 


The ceramic plates are poled tangentially and electrically 
wired in parallel. Thus, it is apparent that, during 
operation, the electric field will either simultaneouly be 
parallel or antiparallel to the poling direction shown 
Figure 3. The ring is wrapped on the outside by epoxy- 
impregnated fiberglass of 0.03 in (0.762 mm) thickness, which 
provides a compressive stress bias. 


A cutaway view is as shown in Figure 4. 
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Figure 4: Cutaway View of Ring Transducer 


The piezoelectric ring and fiberglass wrap is enclosed by 
two encapsulant materials: the inner one is a polyurethane and 
the other is a neoprene. Of particular importance is the 
annular mounting plate at the bottom, which has an inner 
radius less than half the radius of the transducer. Because of 
the narrow size of the center hole of the mounting plate, it 
is expected that the mounting plate will influence the 
acoustic cavity mode of the enclosed water column in the 


Mraadle of the Eransducer. 


Appendix C lists the physical dimensions and constants for 
the materials of the transducer. 

The main operational cnaracteristics of the transducer as 
orovided by the manufacturer are: 

a. Resonance frequencies are 1.0 kHz (cavity) and 2.5 

KHZ (| oanege. 

b. Source Pressure Level, which corresponds to the 

effective pressure on the acoustic axis, at ring resonance 

1S 132 dB ref 1 pPa/V at 1 meter; 

c. Operational Depth, which corresponds to the maximum 

depth where the performance of the transducer 1s not 


compromised - unlimited. 


The measured admittance circles of the ceramic vag 
assembly without and with fiberglass wrapping in air are as 


Shown in Figures 5 and 6, respectively. 
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Figure 5: Measured Admittance Circle of the 
Unwrapped Ceramic Ring Assembly, In 
Sale 
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Figure 6: Measured Admittance Circle of the 
Ring Assembly with Fiberglass 
Wrap, In Air 


The measured transmitting voltage response (TVR) waeeeee 
V drive and the vertical directivity pattern at 2 kHz awemee 


Shown in Figures 7 and 8, respectively. 
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Figure 7: Measured Transmitting Voltage Response 


At the ring resonance frequency of 2750 Hz, the TVR is 132 


dB; it is 127 dB at the cavity resonance frequency of 1250 Hz. 
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The beam pattern is omni-directional in the equitorial 
Dlane and toroidal in a plane passing through the axis of 
symmetry, which in Figure 8 below is the 90° and 270° line. The 
Memiireng plate 1S at 270°. 


The narrow Lobe at 270° 1S probably 
ele 2cO the mounting plate. 


210° 180° 
Directivity Pattern at 2 kHz 


150° 





Figure 8: Vertical Directivity 
Pattern at 2 kHz 
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It. METHODS OF ANALYSIS 


A. INTRODUCTION 

This chapter presents the simple approximate expressions 
for the ring and cavity mode frequencies from G.W. McMahon 
[Ref. 3] and discusses the theoretical background behind the 
finite element code, ATILA, as well as the analytical 
approach, suggested by Hong-zhang Wang [Ref. 2], to obtain the 
approximate electrical admittance equations ORS the 
tangentially polarized piezoelectric thin circular cylindziaeam 


tube in the unloaded eace. 


B. RING AND CAVITY MODE FREQUENCIES [Ref. 3] 

The behavior of the cavity and ring resonance in small 
length-to-radius ratio open tubes of lead zirconate titanate 
ceramic immersed in water was studied by G.W. McMahon. 

When the tube is short i.e., h/a <« nm , G.W. McMahon gives 


the following formulas: 








] 
, ee 2 (1) 
pa° 
eee (2) 
(h + 20a ) 


1 


where w, 1s the ring mode angular frequency, w. is the anguiar 
Frequency of the first cavity mode, ais the inside radius, h 
is the height of the ring, B is the bulk modulus of water, t 
1s the wall thickness of ring and c is the speed of sound in 
Weeer. C, is the velocity of sound in the water column given 
by Equation 5 below and a, an end correction, is given by the 


empirical formula: 


“® = 0.633 -0.106 £22 (2) 


The dimensionless frequency parameter, 9 = w.a/C,, 1S Given by 


a 


oF (pam 
2a 


Gm eot ol 06 O = 5 (4) 


The value of 2 is between 0.33 and 3.3. The velocity of sound 
within the tube, cy, is effectively less than the velocity of 
Bema Cc in open water because of the finite stiffness of the 


meee Wall. It is given by 








The calculated value of the cavity mode frequency for the 


transducer is presented in Chapter V, Section C. 
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C. FINITE ELEMENT ANALYSIS, THE ATILA CODE 

The application of finite element analysis (FEA) to solve 
boundary value problems consists of the transformation of the 
governing differential or integral equation(s) into a2 mama 
nodal matrix equation, the solution of which represents the 
discretized solution of the problem. There are many techniques 
to obtain a finite element formulacion [Reiws 677) 

ATILA is a finite element code developed at 1’Institut 
Supérieur d’Electronique du Nord (ISEN) in France fomaia. 
analysis of underwater transducers. It utilizes the 
Variational formulation of the finite element problem [Ref. 
oO ae oes 

ATILA uses quadratic isoparametric elements. Isoparametric 
means the same polynomial (quadratic in this case) is used to 
interpolate both the geometry and the field variation. ATILA 
has 46 different types of elements, including shell, plate, 
transition, spring, trilaminar, and two- and three-dimensional 
solid elements of various types. With ATILA, it is possible to 
model lossy elastic, piezoelectric, magnetostrictive, magnetic 
and composite materials, fluid, solid-fluid interfaces, 
dampers, and radiation conditions. 

ATILA can perform: (1) static analyses, (2) modal 
analyses, which correspond to a free vibration problem, where 
the eigenfrequencies and eigenmodes are computed, and (3) 
harmonic analyses of radiation or scattering problems, which 


correspond to a forced vibration problem, the excitation being 
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mie Orage sapOlicd across ene Siectrical terminats ®f the 
Seeadsaducer ©r external forces applied to the nodes. 
1. Harmonic Analysis of a Radiating Piezoelectric 

Transducer 

This problem is governed by the equation of motion in 
the elastic and piezoelectric structures, by Poisson’s 
equation in the piezoelectric structures, and by Helmhoitz’s 
Se@ecton in the fluid. Appropriate boundary conditions are 
Semimed, both on the solid-ftiuid interface and on the external 
fluid boundary, which must simulate the appropriate radiation 
Sema 1lion. 


The solid equation of motion is given by [Ref. 12]: 





Ou, do, 
i = i (6) 
: on Ox; 


where p is the solid material density, u is the displacement 
vector, t 1S time, o is the stress tensor, and Ks 1s a 
@eerainate direction. Here i and j can be 1, 2 and 3, and the 
Binstein notation is used, where summation is implied over 


repeated indices in the same term. 


Poisson’s equation is given by [Ref. 12]: 


als: 


where D is the electric displacement vector and x] ewe 
COOrdinate Girection; 2 Canv5eel ovate 
The iinearized,lossiess Helmhoitz Equation for ¥eme 


propagation of sound in fluids is given by [Ref. 14]: 








where V’ is the three-dimensional Laplacian operator, p is the 
acoustic pressure, and t is time. 

In piezoelectric materials the stress tensor and the 
electric displacement vector can be related to the strain 
tensor and electric field vector via material properties by 
the following constitutive equations, which neglect magnetic 


and pyroelectric effects [Ref. 12]: 


where [T] is the stress tensor, [S] is the strain tensor, Eis 
the electric field vector, D is the electric displacement 
vector, [c#] is the constant electric field elastic stiLimese 
tensor, fe] is the piezoelectric tensor, and [e7} is meme 
constant strain dielectric tensor; i, j, kK and 1 can be equal 


EO_1 2 anaes. 


iG 


Ultimately the solution must be expressed in terms of 


ct 


displacements and electric potentials. To chis end the 


—_ 


rollowing two equations from elasticity and electricity, 


—_ 


mesoectively, are used [Rer. 12]: 


sg oe Ou, Ou. . (11) 
d® 
13 = 
| Bs Gy 


where [S] is the strain tensor, u is the displacement vector, 
xX; 18 a coordinate direction, E is the electric field vector 
and ®@ is the electrical potential; i and j can be equal tol, 
2and 3. It must be noted that the quasi-static approximation 
[Ref. 4] 1s used in writing Equation 12 for the low (acoustic) 
frequencies involved here. 

The boundary conditions and prescribed excitations at 
each node can be defined either by a displacement or an 
applied force, an electrical potential or an electrical 
Sharge, Or an acoustic pressure. 

In ATILA, the previous seven equations are transformed 


into the following matrix equation [Ref. 1]: 


[K,,] -w° [M] [K,6] =e) 


ui U F 
[K.o] . [Ko] [0] p = = (13) 
-p?c*w*[L]™ [0]? [H] -w[M] P pc*y 


BLY) 


where: 


Baie: 


vector of the nodal vaiues of the components of 
the displacement field, 

vector of the nodal values of the elect 
votential, 

vector of the nodal values of the pressure field, 
vector of the nodal values of the externally 
applied forces, 

vector of the nodal values of the externally 
applied electricai charges, 

vector of the nodal values of the integrated 
normal derivative of the externally applied 
pressure on the surface boundary S, 

Finite Element stiffness matrix, 

Finite Element pieozoelectric matrix, 

Finite Element dielectric matrix, 

Finite Element consistent mass matrix, 

Finite Element fluid pseudo- stiffness matrix, 
Finite Element consistent fluid pseudo- mass 
matrix, 

Finite Element coupling matrix at the fluid 
structure interface, 

zero matrix, 

angular frequency, 


fluid density, 
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() 


tluid sound speed. 
= means transposed. 


The results of this analysis for each input frequency 
are che complex displacement, rotation, and electrical 
potential fieids at each transducer node, the complex pressure 
field at each fluid node, and the complex electrical impedance 
and admittance. 
2. Modal Analysis of a Piezoelectric Transducer 

This problem is governed by the equations of motion in 
the elastic and piezoelectric structures, and by Poisson’s 
equation in the piezoelectric structures. The matrix equation 
governing this problem is easily obtained from that described 
in the previous section. In a modal analysis there is no fluid 
emer Ehnere are no external forces applied (the natural boundary 


@enditions), so the third row and column of Equation 13 become 


irrelevant, and F is replaced by 0, resulting in 
Ree) — y) [ M] [K,4! | a 2 0 | 1a) 
Bere Kyo) | 2 “2 


where the elements are defined by Equation 13. 

In this equation the resonance condition, which 
@orncesponds to the electrical short-circuit condition, is 
obtained by setting $=0. The anti-resonance condition, which 
corresponds to the electrical open-circuit condition, is 


obtained by setting q=0 [Ref 15]. 
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The results of this analysis are the eigenrrequencies 
and eigenmodes. The maximum number of modes, which must be 


specified By thegusecrawis 0c. 


D. STANDARD NETWORK ANALYSIS [Ref. 2] 
The state equations for piezoceramics as stated in 


Faquation 9 and 10 can generally be expressed as foliows: 


So = Sete ee ee (15) 
Dev= aT 2 Gea as, (16) 
where [T] is the stress tensor, [S] is the strain tensor, E 


1s the electric field vector, Dis the electric displacement 
vector, ({s®] is the constant electric field elastic compliance 
tensor, {d] is the transposed piezoelectric stress tensor, and 
[e7? is the constant stress permittivity tensor; i, j, Kae 
tl can be equal to 1, 2 and 3. The matrix notation has been 
used [Ref. 16]. 

If it is assumed that the ceramic is isotropic in the 


plane transverse to the poling, the piezoelectric equations in 


the cylindrical coordinate system may be written as 
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Sig= SO yy Py AysB. (17) 
Se Sane 


= re , sel oi 
D,=,, T,, + Ay Tyg + Ay, Ti. + €'33-Bp 
Dae Cle ea ee 


where S,,~S33, S)»=S);- If we let the displacement be expressed as 


Pie ,, UU), Enem thewestrain Components are 
du, 
S_= s 
a JU, 7 sea 
= oe 
Cie rE 00 
09 as 
r0é fr (18) 
du, 1: Ou. 
p= t — 
- Boze ae id 
BND 
=" Og 
du. du, 
= = + 
CR tz 
respectively. 


Zale 


The equations of motion of the elastic medium 


cylindrical coordinate system are 














~~ OT). oom ao ete eae 
Oou_= a = nae 
Pc oe al, OZ re 
so ~ WOE nO eg gpI, EASE 
OU= += ~ + 
St Oe oie OZ Lr 
«07... 10 2 ee 
u_= ie _7 Soe 
” Oe at OZ ie 


1. Vibration of Unloaded Tube 


For the tangentially selarizedew nae 


assumptions were made: 


a. the electric field, E 1s approximated as 


(Oe eaes 


b. the unloaded thin-wall tube vibrates with 
axial symmetry 


ae ih 


c. the stress components, T rr 


SAC T'9 , 


and Tg, on the outer, inner surfaces and in 
the whole body of the tube are neglected. 


da fo5 = Olathe ag— en 


£. O1,,/80 = AT,,/00 = JT,,/8 = 0 


ae 


in the 


| Lep 


the foll@waiag 


With the assumptions mentionea, the piezoelectric 


equations may be simplified as 


2 0 
eee ie 5 _ 
liye, = { ve) (1 -u-) : oes (3 7 TA} B, | 20) 
clare at, _ \ : 
ae = 1 soe ean =e . (OS,,7S..— (a, +0d,3} By | 
D, = ys Dy is a, T.. a es 
where Y* = 1/s,,* is the Young’s modulus and o = -s,,*/s,,* is the 
Poisson ratio. 
The Velie een equations OE the thin-walled 
plezoceramic cylindrical tube can also be simplified as 
2 Ee 
pu, = -* 
i 
pu, = 0 (21) 
OT. 
ae ez 
Oe Te 


By combining the piezoelectric equations and the equation 


of motion, the following acoustic wave equation may be 





obtained: 
ee create (22) 
dz : 
where 
ee) 
K? = : a (3) 


Ze 


Sao es : 
£ (a) =:1-¢7} | — =e ; (24) 
ie 8) 


and w, = cC,/a. c. is the acoustic wave speed in the ceramic 


Given by: 


mel 
= Som eae (25) 
p 


Applying the boundary conditions, we obtain the following: 


u. = Aysin(K,z| (26) 


7 6 (dh, + Ody, )+£, (w) (a, +045] 
K,(£,(w) +07] cos 5K, 


J 


a 


a 





If a and b are the inside and outside radii of the img 
and N is the number of silvered electrodes, the total electric 


Charge on the electrodes is given by 


Q = NO, = Mi {, 2 drdz (29) 
-li2 Ja 
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The electric potential between the two near ¢iectrodes is 








V. and from Maxweil's equation, it follows that 
f ie t 
=u an i NVA i \ 
part i—i — 30) 
rij 2m | 


Combining the aforementioned relations, the admittance 


equation are obtained as 








) (31) 
- ae oe 
1-o¢- 
where 
_ Y¥§@ 
He = 31 
's ‘a 
33 
: ; : c32) 
GG lia Fi ao ie ors 





These equations are applied in obtaining the results in 


Chapter V, Section E. 
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IV. TWO-DIMENSIONAL FINITE-ELEMENT MODEL 


As INTRODUCTION 

In order to efficiently modei the free flooded ring, as 
Shown in Figure i, the structure may be considered to be an 
axisymmetric body. Appropriate modifications to the material 
properties must be made to model the structure as 
axisymmetric. A continuing problem in the analysis of these 
plezoelectric devices is the determination of these material 
properties. 

A two-dimensional finite element model was developed using 
8-node isoparametric quadrilateral elements for the 
plezoelectric ceramic, the elastic encapsulant materials, the 
mounting plate, supporting structures, and the fluid. The 
fluid domain extends to 1m from the acoustic center. 6-node 
isoparametric linear elements were used for the solid-fluid 
interface and 3-node linear elements were used for the dipolar 


radiating damper. These were used to terminate the fluid mesh. 
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B. CHARACTERISTICS OF THE MODEL 
1. Material Properties 
a. Unwrapped Piezoeiectric ceramic Ring 

As described before, the transducer has 72 
tangentially poled lead zirconate titanate ceramic staves, 
Glued together with epoxy resin, and arranged in a 0.3668 m 
Gieameler ring. 

In order to simplify the finite element model, an 
equivalent homogeneous piezoelectric ceramic ring with 
Material properties determined by a combination of the 
piezoelectric ceramic and epoxy resin properties was used. 
Because the appropriate combination of properties 1s not 
known, "Smeared" material properties for the homogeneous ring 
were obtained by adjusting the circumferential compliance so 
that the resonance frequencies obtained using the ATILA code 
tolerably matched the in-air experimental measurements of the 
resonance frequencies for the unwrapped ring. 

A limitation of an axisymmetric ATILA model is that 
the tangential poling of the segmented ring cannot be modeled 
directly. It has to be modeled as a radially-poled ring as 


illustrated in Figure 9. 


PAL 


’ 





MODELED AS RADIALLY-POLED RING 





Figure 9: Modeled as a Radially-Poled Ring 


The transformation from the tangentially to radial 
polarization is accomplished by a suitable exchange of the 
elastic, piezoelectric, and dielectric tensors as described in 
the ATILA manual [Ref. 1]. Since the potential degree of 
freedom is a scalar, the change in polarization direction is 
Simply a geometrical transformation. Thus, to model the 
applied potential of “he tangentially poled ring using the 
radially poled ring, the radial potential has to be corrected 
by a factor equal to the ratio of the radial thickness to the 
Stave thickness, d/fe in Figure 9. Similarly, the results 
computed in a harmonic analysis for the radially-poled 
axisymmetric model have to be modified as follows: (1) the 
displacement and pressure magnitudes must be divided by the 


ratio of the circumferential length of a "smeared" 
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piezoelectric element, which is equal to 1/72 of the ring 
circumferential length, to its thickness, and (2) multiply the 
electrical impedance by the square of the same ratio. This 
mero 1552-602 for the transducer under study. 
b. Wrapped Piezoelectric Ceramic Ring 

Normally, the fiberglass belt wrapped around the 
ceramic ring (Figure 3) is modeled as thin shell element. 
However, to simplify the model and to reduce the number of 
degrees of freedom, the 0.762 mm thick fiberglass belt in our 
case was not modeled as a separate element. Instead, the 
material properties of the equivalent homogeneous ceramic ring 
were modified. A sample calculation of the approximate new 
material constants iS given in Appendix D. These approximate 
values were then used for a modal analysis using ATILA and 
were fine tuned by adjusting first the circumferential 
compliance so that the resonance frequency obtained tolerably 
matched the measured values provided by the manufacturer. Then 
the other compliances were adjusted to obtain a close match of 
the coupling coefficient to the measured value. This scheme is 
considered sufficient for the model even though a scheme 
implemented by McMahon and Armstrong [Ref. 17], which fits the 
resonance and antiresonance frequencies and the electrical 
Capacitance taking into consideration the compression bias 
effect of the fiberglass wrapping, might give a more accurate 


result. 


ONS) 


Cc. Encapsulant Materials and Mounting Plate 
The encapsulant materials, namely the polyurethane 
and the neoprene, the reflector ring (high carbon steel), and 
the mounting plate, a fiberglass-epoxy composite, were modeled 
individually as 8-node quadrilateral elastic elements. The 
materlal properties were provided by the manufacturer. 
2. Types of Elements 
The quadratic isoparametric elements listed in Table 
1, which are described in the ATILA User’s Manual [Ref. 1], 


were used: 


Table 1: LYPES Gr ELEMENTS Usa» 


Region Element Geometry 


Piezoelectric ring QUADO8P 8-node quadrilateral 


Encapsulant (Boots, QUADO8E 8-node quadrilateral 
Spacer, Reflector 
Ring) & Mounting 
Plate 


Radiation surface LINEO3R 3-node linear 
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3. Constraints on Mesh Design [Ref. 1] 
Design of the mesh was guided by the following 
constraints: 
a. Aspect ratio 
The aspect ratio of each element should not be 
greater than 3, although 4 is considered an acceptable, though 
less conservative, value. 
b. Internal angles 
The angles between adjacent sides of quadrilaterals 
Should be between 45 and 135 degrees, although 30 degrees and 
150 degrees are considered, respectively, acceptable, though 
less conservative, values. 
c. Element size 
As ATILA utilizes quadratic interpolation 
functions, the size of each element must not be greater than 
one-fourth wavelength at the highest frequency of interest. 
d. Interelement compatibility 
The mesh should be built in such a way that 
adjacent elements have adjoining sides with collocated nodes 
to ensure accurate interpolation at their interfaces. 
e. Radiation boundary elements 
For in-water harmonic analyses (radiation problems) 
the fluid mesh outer limit must be spherical. This is required 
by the radiation damping elements available in the ATILA code. 


ATILA offers damper elements that are attached to the external 


eye 


Surtace of the surrounding fluid domain and are designed to 
absorb various components of the radiated multipolar acoustic 
field to overcome the problem of modeling an infinite fluid 
domain. The ATILA offers the monopole and dipole radiation 
damping elements. The latter includes not only the monopole 
term of the radiated fieid multivolar expansion, but also the 
dipole term. Dipolar damping elements were selected to 
terminate the fluid mesh because, in addition to providing a 
more accurate solution, a smaller mesh can be employed, 
resulting in a lower computational cost. 

A fluid mesh outer limit radius greater than the 
far-field distance is desirable to compute the acoustic source 
pressure level and to compare computed and measured acoustic 
pressure data. The far field criteria for a baffled rigid 


piston-like source is given by [Ref. 15]: 


De 
fe (33) 
and 
eee fe (34) 


where D is the diameter of the ring which is 0.3795 m. 
The minimum wavelength, A, corresponding to the 


maximum frequency of interest is given by: 


A= = 0. iy Sve (35) 


oes 
f 


2 


where c 1s the speed of sound in water, 1500 m/s, and f is the 
mugiest ~-requency Of interest. 4000 Hz, giving Re. > 0.384 m. 

The boundary was placed at a radius R equal to 1m 
from the transducer’s acoustical center, which is more than 
2.5 times the far-field limit of the equivalent piston-like 
aelmctece at the resonance frequency [Ref. 18]. 

f. Boundary Conditions 

Boundary condition data entries enable the user to 
force, for example, clamped, hinged, or simply-supported 
conditions, or to set the master degrees-of-freedom. These 
conditions are generally associated with symmetry planes or 
axes, electroded or pressure release surfaces. 

In a harmonic analysis of a piezoelectric 
structure, the degree-of-freedom, D, for all the nodes 
associated with the applied potential, V, to the plane or 
line, P, containing the node N has to be a master degree of 
freedom and identical. This is realised by using the following 
line in our data file: 

= = 1! V 
for example, 

-l -4 -2 ae 
where the degree of freedom 4 indicates electrical potential 
on a plane 2 (normal to the OY axis) passing through the node 
number 1 with an excitation voltage of 1 V. The negative sign 


in front of the 2 indicates that the electrical potential 
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degree of freedom of all nodes on that plane are identical, 
NG rn ee mY de 

Since our model is a 2-D model, the degree-of- 
freedom in the piane normai to the OZ axis are all deleted 
using the following data line: 

= D ie 
where D is 3 which corresponds to the displacement in the z- 
direction and P is 3 which corresponds to the plane normal to 
the OZ axis. 
4. Final Mesh Designs 

With the above considerations, an in-air 
axisymmetrical mesh about the OX axis was designed for fitting 
the material properties. Only half of the ring cross-section 
needs to be modeled. The bare ceramic staves without the 
fiberglass wrap were divided into five elements to meet the 


aspect ratio requirement of the ATILA code. The in-air model 


is as shown in Figure 10. 
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Figure 10: In-Air Model 


Using the pre-processor mesh generator MOSAIQUE [Ref. 1], 
the in-water mesh without the mounting plate as shown in 
Figure 11 was generated. The transducer mesh consists of 40 
solid elements, 21 interface elements, 166 fluid elements, and 
15 radiating elements, totaling 2884 DOF. This mesh was also 
used to investigate effects of encapsulant material changes. 

The final mesh to be developed was the in-water model with 
the mounting plate. The total mesh is shown in Figure 12 and 
an enlarged view of the solid mesh with the mounting plate 1s 
as shown in Figure 13. There were a total of 87 solid 
elements, 48 interface elements, 406 fluid elements and 35 
radiating elements with a total of 6600 DOF. The input data 


file of this final mesh is given in Appendix A. 
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Figure 11: 
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Figure 12: In-Water Model With Mounting Plate 
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Figure 13: The Enlarged View 
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V. RESULTS 


A. MATERIAL PROPERTIES OF EQUIVALENT HOMOGENEOUS RING 
The material properties of the equivalent homogeneous 
piezoelectric ceramic ring, unwrapped and wrapped, were 
determined by modal analysis using ATILA. This analysis 
corresponds to a free vibration problem in which the 
eigenfrequencies and eigenmodes are computed. 
1. Unwrapped Equivalent Homogeneous Ring 
First, a simple transducer model of only the bare 
ceramic ring without the fiberglass wrap, as shown in Figure 
10, was designed. A modal analysis of this model was conducted 
uSing book values for the piezoelectric properties. Table 2 
lists the calculated resonance and antiresonance frequencies 
and the coupling coefficient for the first ten modes of 
vibration. The coupling coefficient is calculated using the 


relationship [Ref. 19]: 


Meese kK 1s the coupling coefficient, f, 18 the resonance 


Meequency and £. is the antiresonance frequency. 
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TABLE 2: RESULTS OF THE MODAL ANALYSIS OF UNWRAPPED 
PIEZOELECTRIC RING USING BOOK VALUES er 


THE MATERIAL PROPERTIES 


RESONANCE ANTIRESONANCE COUPLING 


HZ ) (HZ) CORFE was 





The ring mode of vibration(mode 2) is the most 
important in operation and is shown in Figure 14. The dashed 
lines in Figure 14 correspond to the rest position, and the 


deformed shape is in solid lines. 
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Figure 14: The Ring Mode of Vibration 


To obtain the smeared properties of the unwrapped ring 
assembly, the circumferential compliance s,, was adjusted so 
that the ring mode resonance frequency obtained for the 
unwrapped ceramic ring in air using the ATILA model matched 
the manufacturer’s measured values. 

Table 3 compares the results obtained using typical 
elastic constant values provided in handbooks and the adjusted 


values. 
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TABLE 3: RESULTS FOR RING MODE OF UNWRAPPED CERAMIC RENG 


USING BOOK VALUES AND ADJUSTED VALUES 









MEASURED USING BOOK ADJUSTEE 
VALUES VALUES 
Resonance OREN Sst | 2491 Hz 235 Saez 
GS 7a (O20 se 
Antiresonance | a2 7 ae : 3287 Hz 3177 Hz 
Soupling 0 ao. lena 






COCBEICI ENE a: 


Note: The number in brackets are the percentage error. 


A comparison of the value for the ring mode frequency 
for the ATILA 2-D in-air model for the unwrapped ceramic using 
the adjusted values for elastic constants and the in-air 
experimental results provided by the manufacturer shows the 
expected match of resonance frequency. The computed coupling 
coefficient, however, was 0.9% higher than the measured value. 
The difference is due to the fact that only the elastic 
constant s,; waS adjusted. (Several ATILA modal analysis runs 
were made with adjustments made to the elastic constants sj, 
and S, just to see the effect on the ring resonance frequency 


and coupling coefficient. The resonance frequency changed only 
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peeSsOuevesy cMail.;, Tt was ceit chat a 0.9% diftterence is 
Heeeptable TOr our purposes. For our results, the value of Sy, 
was increasea by 12 percent i.¢e., more compliance) from the 
Beem Value of 0.1539E-10 to 0.1724E-10. This 1s consistent 
with our expectation that a ring made of glued segments is 
fiemee COMDliant than a solid piezoelectric ceramic ring. This 
value of s,;, was used for calculating the first-guess smeared 
material properties for the fiberglass wrapped ceramic ring as 
described below. 
2. Wrapped Equivalent Homogeneous Ring 

Smeared material properties for the homogeneous 
equivalent wrapped ceramic ring assembly were obtained by a 
four-step process. First, the coarse values for the elastic 
constants were calculated. A sample calculation is given in 
Appendix D. The circumferential compliance s;, is then adjusted 
so that the calculated ring resonance frequency matched the 
measured frequency. The other elastic constants were then 
adjusted by a common factor such that the computed coupling 
coefficient matched that of the measured value. The results 
for the ring resonance frequency and coupling coefficient are 


given in Table 4: 


43 


TABLE 4: RESULTS OF THE WRAPPED EQUIVALENT CERAMIC RING 


OBTAINED VALUES 


. : . 
MEASURED USING VALUES | ADJUSTED 
| 
! 


RARLIER 


RESONANCE 2462 Hz 


FREQUENCY 


COUPLING 


COEFFICIENT 





Note: The number in brackets are the percentage error. 


Table 5 lists the values of the elastic constants 
calculated as described in Appendix D, and the final adjusted 
values. These latter values were then used for in-water 


harmonic analyses: 
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Pesiobe os MoUASTiC CONSTANTS USED FOR IN-WATER HARMONIC 


ANALYSES 


COMPLIANCE COMPUTED AS IN | ADJUSTED .m°/N) 


| 


CONSTANTS APPENDIX D (m*/N) 


Sis= Ss 3.9e-11 


S66 Se@25eE7L7 See] i 





The other material constants, such as the dielectric 
constants and the piezoelectric constants, were not adjusted. 
The effective density of the homogeneous ring incorporating 
the fiberglass wrap was obtained by dividing the total mass of 
the ceramic and fiberglass by the volume of the homogeneous 


€-ramic ring. 
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Bs IN-WATER HARMONIC ANALYSES 

A harmonic analysis corresponds to a forced vibration 
problem, the excitation being the voitage applied across the 
electrical terminals of the transducer. The applied voltage in 
our case 1S 1 Vrms. Internal material losses for the 
encapsulant materiais were included in this model (but not for 
the ceramic ring assembly). These were obtained from the loss 
tangents for the respective materials. Values of 0.4 and 0.7 
were used for the loss tangent for polyurethane and neoprene 
respectively [Ref. 20]. 

Ll. Mounting Plate Not Included 

Using the material property values obtained from in- 
air modal analyses, described earlier, an in-water model 
without the mounting plate was first analyzed. The mesh is 
shown in Figure 11. 

The plot in Figure 15 depicts the transmitting voltage 
response curve obtained by ATILA for this model along with the 
corresponding manufacturer’s data for the actual transducer 
(with the mounting plate). 

The computed transmitting voltage response was 
generally of the same form except for the lower frequeuay 
region between 1000 to 1750 Hz. The discrepancy in this region 
is probably due to not including the mounting plate in this 
model. The mounting plate is expected to change the 


transmitting voltage response at the cavity resonance as well 
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aS Une Cavity resonance frequency itself. The TVR at the ring 
Be-OUance Exequency was Calculated to be about z dB lower than 


that measured at the same frequency. 
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Figure 15: Comparison Of The Measured Transmitting Voltage 
Response With That Computed For In-water Model 
without Mounting Plate 
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2. Mounting Plate Included 
An in-water ATILA model including the mounting plate 
was then developed. The mesh are shown in Figures 12 and 13. 
A comparison of the calculated transmitting voltage response 


with the measured response is shown in Figure 16 below: 


140 , 








— Measured 


dB re 1uPa/V at 1 Meter 






—+ Computed 


100 7 7 
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Figure 16: Comparison Of Measured TVR With That Calculated 
For In-water Model with Mounting Plate 
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The agreement between the calculated and measured TVR 
iS about the same as the previous model. It is creatly 
improved in the cavity mode region and the ring resonance 
frequency calculated is exactly the same as the measured 
values. However, there is a deep dip at 1500 Hz. If this dip 
can be made shallower, then the overall agreement would be 
excellent. It is hoped that the results will be better when 
accurate material property values are available for the 
polyurethane and neoprene encapSuiants. The variation in TVR 
with encapsulant materiai properties is explored in the 


Section D. 


C. COMPARISON OF CAVITY MODE FREQUENCY [Ref. 3] 

The results of the previous section for the cavity mode 
frequency can be compared with the approximate analytical 
formulae given by G.W. McMahon and present previously in 
Saapeer Lil. 

Using relationships given by G.W. McMahon, the first 
cavity mode can be calculated from Equation 2. The speed of 
sound in the water column, cq is first computed using Equation 
5 presented in Chapter III. With B equal to 2.18 x 10° N/m and 
fogeeo 5.975 x 10° N/m’, c, is found to be 846.06 m/s. 9 is 
found to be 1.075 and @a is 0.519. With these values, the 
cavity resonance frequency, w,, iS 7320.7 rad/s or f, is 1165 
Hz. The relative error between this value and the value 


obtained from the modeling (1250 Hz) is 7%. 
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D. EFFECT OF ENCAPSULANT MATERIAL PROPERTIES CHANGE 

The ATILA model, without che mounting plate, was used to 
explore the effect of changing the material properties of the 
encapsulant materials. namely cthe polyurethane around the 
ceramic and the neoprene boots. 

Figures 17 and 18 show the results of the TVR computed for 
different values of the Young’s modulus for neoprene and 
polyurethane, respectively. The curves are plotted over the 
range of values oof Young’s modulus quoted in the 


Manltacturer Ss Viterarune: 
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Figure 17: Effects of Varying Properties of Polyurethane 
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Pigure 18: Effect of Varying Properties of Neoprene 


From Figure 17, it can be seen that the best agreement 
with the measured result is for the curve with Y = 4.14 x 10° 
N/m and in Figure 18, it was the curve with Y = 1 x 10’ N/m’. 
These two values were used for the finite element in-water 
harmonic analyses presented in the previous sections. 

It is observed from the two figures that the value of 
Young’s modulus for polyurethane has much greater influence on 
the TVR than that for neoprene. No noticeable change in ring 
resonance frequency was observed in any case. The transmitting 


voltage response was found to increase with a decrease in the 


ayib 


Young’s modulus of polyurethane. It was, Nowever, the reverse 


in the case of the neoprene. 


E. THE ANALYTICAL MODEL 
1. The Unloaded Case 
The expressions developed based on work by Hong- zhang 
Wang [Ref. 2] were implemented using the MATLAB program. The 
program iisting is provided in Appendix B. The effective 
elastic constant values obtained for the unencapsulated, 
fiberglass wrapped ceramic ring assembly of the finite element 
model were used for the analytical model. The admittance 
circle of the unencapsulated fiberglass wrapped ring is 


graphically compared to the measured values in Figure 19. 
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Figure 19: Comparison of Analytical with Measured Results of 
Wrapped Ring In Air 
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The resonance frequency obtained using the analytical 
Metnoad (2420 Hz) was approximately 2 percent lower chan the 
measured (2462 Hz) and the value obtained using che finite 
eeement code (2465 Hz). 

Wang's approach is quicker and easier to use compared 
to finite element modeling as only the circumferential 
compliance s,, and the Poisson ratio are required for the 
modeling. It is especially useful when oniy measured 
admittance circle and handbook material property values are 
available. The expressions developed, however, do not take 
into consideration the encapsulant materials. 

Also, Wang’s approach could not be applied to the 
beezoclectric thin ring transducer in water, as it was not 
possible to obtain expressions for the radiation impedances 
presented to the inner and outer surfaces of the ring 


transducer. 
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VI. CONCLUSIONS 


A successful two-dimensional finite element model of a 
free-flooded, tangentially polarized plezoelectric ring Sema 
transducer was built. The model includes internal losses of 
the encapsulant materials. Although the model includes many 
necessary simpiifications to handle the problem in the 
available MICROVAX VMS system, they were designed in such a 
way that the transmitting voltage response results obtained 
from the in-water harmonic analyses differs by only 2 dB from 
the measured value at the ring resonance frequency. In the 
region of the cavity mode, the model including the mounting 
plates has better agreement except for a dip at 1500 Hz. It 
was able to accurately predict the ring resonance frequency. 

The model built (without the mounting plate) was used to 
investigate the effect on the acoustical properties of the 
transducer due to a change in the properties of the 
encapsulant materials used in manufacturing. It was observed 
that the Young’s modulus of the polyurethane has a significant 


influence on the computed transmitting voltage response. 
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For the anaiytical modei, only the case of the transducer 
in alr was developed. The resonance frequency was found to be 


only 2 percent lower than the measured value. 


a1) 


VII. SUGGESTIONS FOR FUTURE WORK 


Suggestions [or futume work sanclude. 

1. Obtain samples of the actual encapsulant materiais 
used for testing to determine the actuai material 
properties. With these vaiues, the model can then be 
refined. A harmonic analysis should then be performed 
and the transmitting voitage response results obtained 
compared with measured values. 

2. With known values of the material properties, the 
effects of changes in encapsulant materials should 
again be investigated using the in-water model with 
mounting plate. 

3. Finally, it would be very helpful and interesting 
to obtain an analytical expression for the radiation 
impedance of a free-flooded ring transducer. With this 
and the results of Hong-zhang Wang [Ref. 2], a more 
complete analytical model for the transducer in water 
can be produced. The results obtained with this model 
should then be compared to those obtained using ATILA 
and the measured values. 

4. Study of array element interaction effects should 
be made using these models to develop techniques for 


array performance modeling. 
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APPENDIX A 


INPUT DATA FILE FOR HARMONIC FINITE ELEMENT MODEL 


% In-Water Harmonic Analysis of a free-flooded piezoelectric ring 
ei sonar transducer 

- Prepared By MAJ Tiong Bemg Tay, .14° Feb 1993 

* 


——_—— ee ee ee ee eee eee ee ee eee ee ee eee ee ie ie ee eee ee ee ee ee ee eee eee eee 


* Navy Type I ceramic ig 
ei Encapsulant Material: Polyurethane & Neoprene a 
* Tangentially Polarized but modeled as radially-poled * 
* Radius of dipolar damper = 1m e 
Woeomao op a SS SS SS SS a a SS St ea ee ee * 


——_- — cc wc ww ew ee ee ee ee ee eee ee ee ee ee ee ee ee ee ee ee ee 


RADIATION DIPOLAR 
ANALYSIS HARMONIC 
SKYLINE COMPLEX 
PeeciolON DOUBLE 


CLASS AXISYMME 
EerpDC = 7 
NLOAD = 40 * The number of loading cases must be less than or equal 
* to 100. The minimum is the number of frequencies 
* multiplied by 2 because we are analyzing a radiating 
* transducer 
FREQUENCY 
0.225E+04 0.250E+04 0.275E+04 0.300E+04 0.325E+04 0.35E4 & 
O.37SE4 0.4E4 
GEOMETRY POLARIZA CARTESIA 
a 
0.900E+02 0.000E+00 0.180E+03 
GEOMETRY 
Bs 
0.762E-03 0.000E+00 
3 
Q0O.1L00E+01 
MATERIAL 
PZT4TA 
peeUOE+O0 O0.Q000E+00 0.700E+04 0.000E+00 0.000E+00 0.000E+00 & 
QO.134E-10 -0.443E-11 -0.405E-11 0.000E+00 0.000E+00 0.000E+00 & 
-0.443E-11 0.170E-10 -0.443E-11 0.000E+00 0O.0O00E+00 0.000E+00 & 
-0.405E-11 -0.443E-11 0.134E-10 0.000E+00 0O.000E+00 0.000E+00 & 
Q0.Q000E+00 0.000E+00 0O.000E+00 0O.390E-10 0.000E+00 0O.000E+00 & 
0.Q000E+00 0.000E+00 0.000E+00 0O.000E+00 0.357E-10 0.000E+00 & 
O0.Q000E+00 0.000E+00 0O.000E+00 0O.000E+00 0O.000E+00 0.390E-10 & 
O0.Q000E+00 0O.000E+00 0O.000E+00 0O.000E+00 0O.S525E-09 0.000E+00 & 
0.000E+00 0.000E+00 O.000E+00 0O.525E-09 0O.O000E+00 O.000E+00 & 
-0.135E-09 0.300E-09 -0.13SE-09 0O.000E+00 O.000E+00 0O.000E+00 & 
0.624E-08 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 & 
O.Q00E+00 0.624E-08 0O.Q00E+00 0O.000E+00 0O.O000E+00 0O.O000E+00 & 
0.000E+00 0.000E+00 0O.S70E-08 0.000E+00 0.000E+00 0.000E+00 
FGLR 
O.587E+11 0.280E+00 0.254E+04 
jealisysh®) 
O.1LQ00E+09 0.490E+00 0.110E+04 0.000E+00 0O.000E+00 0.400E+08 
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NEOPRENE 
0.100E+08 
HRS 
O -1SSE+i2 
KENTRAN 
O-135E+05 


oa 
O 
8. 
es) 
ep) 


WON AU PWN +b 


+ + + + + £ + ££ F F F FF OF F OF ~ F OF OF OF ~ F OF OF « F~« FF: F * OF « OF «* OF «* OF«O F~« F«*OF:~* OF :O OF «OF * FOF :~* F «OF OF * F~*OF * OF OF: F FO KF OF: OF :* OF «OF ~ OF ~ OF 
NO 
oo 

+ + + © + & FF + FF F + F F FF F FF F OF OF «* OF OF « FF * OF « OF «* F OF ~ OF « OF * OF «* OF« OF «OF « OF «OF * OF «OF ~« OF ~* F * OF * OF «* OF« OF ~*OF~«O OF * OF ~* OF *OF~ OF OF * F * OF ~* OF *OF 


e 


See ene GS) CS) OOO OO OO OC OC Oe ae 


0.499E+00 


0.333E+00 


0.330E+00 


(ebog0b-O1 
-s6se0n-01 
-cGo350R 400 
2755658 0m 
oe Gomi 
.64770E-01 
-64770E-01 
-64770E-0O1 
2939/7 S5e20 1 
SSG) (Siga 0d 
45 250m 01 
-43180E-01 
-43180E-01 
os oom-O 1 
so2325e 00 
a2 So Oe oa: 
yolS9 0 = On: 
s2lSf0R- OL 
07958 -021 
ul 7 or] Onl 


OO000E+00 


- OO000E+00 
-0O0000E+00 
SLOTS Se oF 
- 107958 02 
-215908 021 
o2 050 080 
-2LsceE- 08 
sieoe St Ow 
S25 55801 
-43180E-01 
-43180E-01 
.43180E-01 
1530 2 on Ou 


539 755-02 


.64770E-01 
-64770E-0O1 


64770E-01 
ISSeSe-OL 


ts ys) 5 sys)— 1 0)at 
.86360E-01 
7o6so00B- 0 
.86360E-01 
Ae) 3 
-0.15875E+00 
=O .15375E+00 
-0.15240E+00 
-0.15240E+00 
-0.14605E+00 
-0.14605E+00 
-0.14605E+00 
-0.13350E+00 
-0.13350E+00 


15875E+00 


0.148E+04 


0.770E+04 


O0.170E+04 


.18339E+00 
-18656E+00 
.18974E+00 
.18339E+00 
.18974E+00 
.18339E+00 
.18656E+00 
.18974E+00 
.18339E+00 
.18974E+00 
.18339E+00 
.18656E+00 
.18974E+00 
.18339E+00 
.18974E+00 
.18339E+00 
.18656E+00 
.18974E+00 
.18339E+00 
.18974E+00 
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0.16942E+00 
0.17374E+00 
0.16510E+00 
0.17374E+00 
0.16510E+00 
0.16942E+00 
0.17374E+00 
0.16510E+00 
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a) ee Ca eee Con) () 6) eS) Se) OO) Oe OO Oe e oOo OO OO Ooo Oo 0 0 OO 2 
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0.00000E+00 
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0.Q00000E+00 
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APPENDIX B 


MATLAB PROGRAM FOR ANALYTICAL MODEL OF UNLOADED TRANSDUCER 


% 

*s This MATLAB program was written by MAJ Tay Tiong Beng 
’ on 1 Feb 1993 to implement the formula developed 

* by Hong-zhang Wang [Ref. 2] for unloaded tangentially 
% polarized ring transducer 

= 

6. Jn 1 etalizacionu 

clear; 

Gig. 

'del res 


'del resl.met 
'del res2.met 
'del res3.met 
diary res 


* Material Constants 


j=sqrt(-1); 
rhoc=7002; 
s11=0.1702e-10*(1-0.015*j3); % Losses is 1.5% of the real part. The real part 
% is the published value of the elastic constant 
Sibz=-0.405e-015, 

nu=- (s12) /(s11) % Poisson Ratio 

s66=2* (s11-s12) 
e33=0.57e-8 
d31=-0.135e-9 
a33=0.3e-9 
Ye=1/sl1l; 
F=3000; 


oe 


Frequency 

© Fluid Constants 

Calr=343; 

cwater=1500; 

CCer=sGne (ue, Hnee)y, 

* Dimensions 

Inside radius 


outside radius of ceramic stave 
height 


a=(14.44/2) *0.0254 % 
b=(14.94/2) *0.0254 % 
l=658*07 0254 z 


% Capacitance and Admittance 


N=72 ? 

Co=N*2*e33*1/ (2*pi) *log (b/a) 

k31s=Yc* (d31°2) /e33 

t= (d33/d31) 

Cs=Co* (1- (k31s/ (1-nu%*2) ) * (1+2*nu*t+t*2) ) 


o6 


Peeeaamittance Circle 


J=saqrt(-1) ; 
wr=ccer/a * Radial Res 
meaewit/ (2*p1) 
nN=0; 
femoral: S:F 
fa—T1+1 ; 
Wet ft: 
WW=W/WE; 
WW (1,n) =WW; 
Jwcs=]*w*Cs; 
fa=(l-nu'2)*(w/wr)*2 - 1; 
ew CO*K31s/fa; 
ka2=(w/ccer) “2*fa/((w/wr) *2-1) ; 
Rael) =kaZ2 ; 
ka=sqrt (ka2) ; 
hal nn) =ka; 
B= (fa* (l+nu*t) +nu~ (nut+t)) *2/((1-nu*2) *2* ((w/wr) “2-1)); 
B=(nutt) 2/(1-nu’ 2) ; 
meee (tan(0.5*ka*)) /(0.5*ka*1l)) -D)+jwes; 
mm / {4 *wr*Co) ; 
dam a gg 
Za) =Z 3 
elem) =) > 
beers 1) =TT; 
end; 
q=abs (Y) 
G=real (Y) 
S=imag (Y) 
le sb 5 


axis (’square’) ; 


meet (G,>s) ;grid; 

title (’ADMITTANCE CIRCLE OF A TANGENTIALLY POLARIZED UNLOADED TRANSDUCER’ ) 
xlabel (’Conductance, G mhos’) ; 

ylabel(’Susceptance, S_ mhos’); 

meta resl; 


wteeif .G); grid; 

title (’CONDUCTANCE AS A FUNCTION OF FREQUENCY’) ; 
xlabel(’Frequency, Hz’); ylabel(’Conductance, G mhos’); pause; 
meta res2; 


Bemeeic’,S); grid; 

title (’SUSCEPTANCE AS A FUNCTION OF FREQUENCY’) ; 
xlabel(’Frequency, Hz’); ylabel(’Susceptance, S mhos’); pause; 
meta res3; 


diary off 
end; 
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APPENDIX C 


TYPICAL VALUES FOR MATERIALS CONSTANTS 


TABLE 1: PHYSICAL DIMENSIONS AND MATERIAL CONSTANTS 





The Transducer 


Inside Radius ee Cees sys) ical 
Outside Radius “20.2072 4am 
Height 2702292 aan 


Motor Element: Navy Type 1 Ceramic (Hexagonal Class 6mm) 


Density : 7550 kg/m 
Inside Radius [02 8238 3am 
Outside Radius . 0 138973 som 
Elastic Constants (m’/N) Piezoelectric Constants 
g tensor (m’/c) 
Sos eon ad censor m/e 
E — -_ - 
Sj, = -5. 318-12 pee 
E —w. = = 
s— = -4.05E-12 SS it ee 
E = _ 
SHE: = tee oimecia a, = 40.5E-3 
Su =) oe OH aa 
Scgy = 3.25E-11 Giemsa e rm main a7 
Dielectric Constants d,, = 300E-12 
3, = -135E-12 


€3,5 = 1.149E-8 
€;,) = 1.2963E-8 
€,;, = O0.57E-8 
0.624E-8 


dhs = 525E-12 


Mm 
—_ 
= 

tt 
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TABLE 1: CONTINUED 


Fiberglass Wrap 


Thickness 02762 mm 

Density : 2540 kg/m 
Longitudinal Modulus 8.5E6 psi (5.86E10 Pa) 
Transverse Modulus 2.9E6 psi (2.00E10 Pa) 
Poisson Ratio Ons 


Encapsulant Material 


Polyurethane 


Density :0.039 pci ( =1100 kg/m} 
Young’s Moduius : 600 psi (4.14 MPa) 
Poisson Ratio sO 249 

Neoprene 

Density : 0.045 pci ( #1250 kg/m’) 


Young’s Moduius : 100-3000 psi 
(0.69-20.69 MPA) 
Poisson Ratio : 0.49 


Reflector Ring: HRS Carbon Steel 


Density 7700 kay im 
Young’s Modulus : 19.5E10 Pa 
Poisson Ratio SOs 


Mounting Plate: Fiberglass-epoxy composite 


Density : 1700 Kg/m 
Young Modulus : 2000 psi (13.79 MPa) 
Poisson Ratio : 0.34 


ae 


APPENDIX D 


SAMPLE CALCULATIONS FOR THE EFFECTIVE 
ELASTIC CONSTANTS OF THE FIBERGLASS WRAPPED CERAMIC 
RING 


The effective elastic constants of the ring are estimated as 
follows. Consider a section of the fiberglass wrapped ring as two 
rectangular volume elements cemented together as shown in Figure 1 
below. The purpose is to calculate the elastic constants of an 
equivalent ceramic material which will give the same elastic 
properties as the fiberglass wrapped ceramic ring. The equivalent 
ceramic material has the same length, 1, and height, h, as the 
fiberglass and ceramic, and a width equal to the sum of the widths 


of the fiberglass and ceramic. 





Figure 1 


100 


The thickness of the fiberglass is t and that of the ceramic is w, 
where t and w are 0.000762m and 0.0064m, respectively. 

Semeider first an applied force, F,, in the circumferential 
direction only, as shown. Assuming that the resultant stress, T,, 


meme cOorm then the strain, S,, in each material is given by: 


S, = =a = sce, (1) 
For the ceramic, 

(Azz) . A,( S13). T; (2) 
And for the fiberglass, 

(An) , = A, (S813) , He (3) 


There is no reason to assume (Ah) equals to (Ah)... We 


g 
mg@ererore define (4h) .-- for the equivalent material as the average 


Ah across the composite material, neglecting boundary effects: 


(Ah) .w+(Ah),t = (Ah),, (w+t) (4) 


g eff 


Then the effective strain in the axial direction, i.e., in the "1" 
direction is given by: 


yp (5) 


@eees taking h. = h. = h. 


LoL 


Thus the effective CompiManes ic ets S 6) oes 


So i lee ie 
\513/ , (S13), es 
ee eG 


\ S,3) et 7 
By considering forces in the other directions, im ama 
Same result is obtained, with the appropriate identification of 


indices and element dimensions. For exampie, to compute s,, we will 


consider an applied force in the "2" direction, 1.e€., theme 
direction, and calculate the strain in the "1", Gi 
digger lon. 


Using values listed in Appendix C for the fiberglass and 


ceramic, the following were obtained: 
CONSTANTS COMPUTED VALUES (m’/N) 


to a Joe ee Oe 


= 2 eae Oa 
2593 woe Ose 


= 33 a0 
it, FHOw) oe ager 
27) Be ga) 


a 5 Gao er O! 
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